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'^ ■ Abstract 

The most general two Higgs doublet model contains new sources of flavour violation 

(N ■ 

^ ■ that are usually in conflict with the experimental constraints. One possibility to suppress 

o 



the exotic contribution to the flavour changing neutral currents consists on imposing the 
r;; ; alignment of the Yukawa couphngs. This condition presumably holds at a high-energy 

• I scale and is spoiled by the radiative corrections. We compute in this letter the size of the 

O ' radiatively induced flavour violating Higgs couplings at the electroweak scale. These 

O ! 

also yield the absolute lower bound on the size of the exotic contributions to the flavour 

changing neutral currents in any two Higgs doublet model, barring cancellations and the 

^ . existence of discrete symmetries. We show that these contributions ai^e well below the 

^ . experimental bounds in large regions of the parameter space. 



1 Introduction 

The LHC finally set out to find the Higgs and in anticipation of its results it is important 
to consider all viable scenarios for the Higgs sector in order to be able to interpret the data 
once they arrive. The simplest extensions of the SM Higgs sector are so-called two Higgs 
doublet models (2HDM) where a second Higgs doublet is added to the SM particle content. 
In addition to the three Goldstone bosons needed to generate the gauge boson masses, there 
are then five physical Higgs particles: two neutral scalars h and H, one neutral pseudo-scalar 
A and two charged scalars H- . 

In the most general 2HDM there are two Yukawa matrices per fermion type (up-type 
quarks, down-type quarks, leptons) which cannot be simultaneously diagonalized. Without 
any further protection, this leads to unacceptably large flavour changing neutral currents 
(FCNCs). Therefore in the literature mostly 2HDMs with a discrete symmetry have been 
considered, ensuring that each right-handed fermion field can couple to no more than one 
Higgs doublet. There are several choices of the discrete symmetry. The ones most often 
studied are referred to as type I (only one of the Higgs doublets couples to the fermions) or 
type II (down-type quarks and leptons couple to one Higgs doublet, up-type quarks to the 
other) 2HDM. FCNCs are completely absent at tree level in such models. 

A generic way to suppress the FCNCs consists in imposing the hypothesis of Minimal 
Flavour Violation to the flavour symmetry breaking parameters [1]. This hypothesis is im- 
plemented in the aligned two-Higgs Doublet model recently discussed in [2]: If the two 
Yukawa couplings for each fermion type are aligned, they are simultaneously diagonalizable 
and FCNCs are absent at tree level. This ansatz is more general than discrete symmetries in 
two ways: Firstly, it contains type I and type II and all other 2HDMs with discrete symme- 
tries as particular cases. Secondly, it allows to have additional sources of CP violation in the 
Yukawa sector, in contrast to type I and II 2HDM. 

The alignment condition presumably holds at a high-energy scale and, in general, will 
be spoiled by quantum corrections. In this letter we analyse the size of these corrections in 
order to determine the viability of this scenario. Furthermore, in the absence of cancellations 
and in the absence of ad hoc discrete symmetries, this scenario yields the lower bound on 
the exotic contributions to the FCNCs in any 2HDM. 

The paper is organized as follows: In sec. 2 we recapitulate the 2HDM and Yukawa 
alignment. In sec. 3 we present an approximate solution to the renormalization group equa- 
tions (RGEs) in a 2HDM with Yukawa alignment. Therefrom we derive expressions for 
the flavour violating neutral Higgs couplings in the quark sector in sec. 4. Bounds on the 
parameters entering in these couplings are then derived from experimental constraints on 
meson- antimeson mixing and the leptonic B decay in sec. 5. Finally we conclude in sec. 6. 



2 The 2HDM with Yukawa ahgnment 

The general 2HDM we consider has the same fermion content and SU(3)c x SU(2)l x U(l)y 
gauge symmetry as the Standard Model. The Higgs sector consists of two scalar SU(2) dou- 



2' 

and the Yukawa part of the Lagrangian therefore reads: 



blets 01 and 02 with weak hypercharge Y = j. Both Higgs doublets couple to all fermions^ 



+iY'u%q'LiURjh + {Yf\q',4'^jcf>2 + (yf )),yl,4/2 + h.c. 

where i,j are flavour indices and ^a = iT24>*a- The neutral components of the two Higgs 
doublets acquire vacuum expectation values (vevs) during electroweak symmetry breaking 
(EWSB) which in general can be complex. While one phase can be rotated away, the phase 
difference is physical. Nevertheless we can choose to work in a basis where both vevs are 
positive, shifting the phase to the potential and the Yukawas: 



<0«) = ^ 

V2 





V Va 



(2) 



In order to minimize the FCNC effects without imposing ad hoc discrete symmetries we 
postulate that at a high energy cut-off scale. A, the Yukawa couplings of the same fermion 
type are aligned [2]. We parameterize this condition as: 

Yi'\A) = cos iA„F„, yf\A) =siniA„F„, (3) 

Y^/\A) = cos ifi.Y,, Yf\A) = sin 0,7^, (4) 

Y^'\A) = cos iA,y„ yf^(A) = sin 0,7,; (5) 

The type I 2HDM is contained in this parameterization as the special case 0^ = 0£/ = 0, = 
and type II as 0„ = 0, 0^ = 0, = |. 



3 Radiative corrections to the ahgned Yukawa coupHngs 

The renormalization group equations (RGBs) for a general 2HDM with Yukawa alignment 
have been derived in [3]. We reproduce them in our notation in appendix A. The radiative 
corrections introduce a misalignment of the Yukawa couplings at low energy. To see whether 
this leads to unacceptably large FCNCs, we solved the RGBs numerically and analytically 
using the so-called "leading log approximation" which estimates the down-type quark cou- 



'This is sometimes called a type III 2HDM. 



plings at the electroweak scale as: 



Ak)i 



Ak) 



(mz\ 



r-{mz) . i^-(A) + ^^^^.,(A)log^-^j 



(6) 



and similarly for the Yukawa matrices of the up-type quarks and leptons. Inserting the yS- 
function (49), the coupling at the EW scale takes the form: 



rik). 



r,«l 



Y'^\mz) « k'^'Y, + efy^jlY, + 6fY,YlY,, 



(7) 



where the coefficients /r^ , 6^ ^ and ^^^^ can be found in appendix B, as well as the correspond- 
ing formulae for up-type quarks and leptons. 



4 Flavour violating neutral Higgs couplings 

To derive the low energy Lagrangian it is convenient to rotate the Higgs fields to a basis 
where only one of the two doublets, say Oi, gets a vev: 






2 ; 



' cosjS sinjS 
sinyS -cos/3 ^ 



02 , 



(8) 



Here tan/3 = ^ is the ratio of the vacuum expectation values in the original basis. In the new 
basis the Lagrangian can be written in the following form (quark sector): 



XYukawa = |^L (M„Oi + T^Oj) u',^ + ^^ (M^Oj + T^/Oj) ^ + h.c.) , 



(9) 



where v^ = v\ + v\ = (246 GeV)^ and the couplings Mu4 and Yu4 are evaluated at the scale 
ruz. Their expression in terms of the original couplings in the basis {^1,02} is: 



Md,u{mz) = -^ (cos/3 Y^limz) + sinyS Y^^{mz)) , 
Td,u{mz) = ^(- sinyS Y^^l{mz) + cos/3 Y^imz)) 



(10) 
(11) 



In order to rewrite the Lagrangian eq. (9) in terms of the mass eigenstates, we first express 
the Higgs doublets Oi , O2 in terms of the physical Higgs fields h,H,A, H- and the Goldstone 
bosons G°,G-: 

^ -^{v + cos(a - p)H - s\n{a - /3)h + ?G") J ' 



O, = 



(12) 



o. 



(13) 



■^(sm(a - /3)H + cos(a - /3)h + iA) 

for a CP conserving Higgs potential and where a is the mixing angle of the mass eigenstates. 
Note that in a general 2HDM the ratio of the expectation values tan^S has no well defined 
meaning. The basis of the Higgs fields can be freely chosen and we could just have started 
in the basis Oi, O2 instead of ^1, 02 (thus setting /3 = 0). The only relevant mixing angle is 
therefore a - /3. This is in contrast to 2HDM type I and II where there is a clear distinction 
of the two Higgs doublets by the way they couple to the fermions. The ratio of the two vevs 
then gets a real, physical meaning. 

Finally, we perform unitary transformations in flavour space of the quark fields 

"l = yu(.fnz) Ul, Mr = Vj^imz) Ur, (14) 

4 = y,^(mz) dL, 4 = ydimz) dR, (15) 

in order to diagonalize the quark mass matrices: M„ = V^M'^'^^'VuK M^ = Vj^Mf^^'V^. In 
this new basis, where the Higgs and quark mass matrices are all diagonal, TiXmz), ^di^z) 
are not diagonal and thus give rise to the following flavour violating neutral Higgs couplings: 



X 3 Ml Am [cos(a - /3)/i + sin(a - /?)// - zA] ur 
+ diAd [cos(a -/3)h + sin(Q' -fi)H + iA] dR, 



(16) 



where: 



K = ^Vl;\mz)T„(mzW!!(mz), (17) 

Arf = -Vf(mz)r,(mz)y^(mz). (18) 

It is possible to calculate approximate expressions for A„, A^ noting that: 

A„ = V„'^r„M-iy„^)(y„^+M„y„«), (19) 

and analogously for A^. Substituting eqs. (7), (10) and (11) and keeping the lower order 
terms in e,,, 5u we find that the ofi'-diagonal couplings read: 

^"''- = EMu, (20) 
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(b) Numerical result for 2.5Ei 
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(a) Analytical approximation of E^ 

Figure 1: Contour plots ofE^forK - lO'^ GeV. The left figure corresponds to the analytic formula, 
eq. (25). Solid/dashed/dotted lines correspond to the absolute values of 1/0.3/0.1, blue lines corre- 
spond to negative values, green lines to positive ones. The right figure shows 2.5Aj_23/2rf,23 where 
Arf,23 hcis been obtained by numerically solving the RGBs. The rescaling was done in order to make 
the comparison to the analytical result easier. 



where, assuming real i/r„, \J/^f. 
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(22) 
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(24) 
(25) 



Thus, the ofF-diagonal elements of the flavour violating Higgs couplings A„ ,/ can be factor- 
ized in two parts: Q^^ are determined by the experimental values of the entries of the CKM 
matrix and the quark masses, whereas Eu4 depend on the unknown details of the 2HDM and 
the scale A at which the alignment condition is imposed. It is apparent from (23) and (25) 
that Eu4 depend just on two parameters: jS - i/r„ and /3 - if/d- Moreover, for i/r„ = ifj^ and 
ipu = iffd ± 7r/2 the flavour violating Higgs couplings vanish, since £„ = E^ = 0. This choice 
includes as special cases the type I (i/r„ = if/^ = 0) and type II (i/^„ = 0, ij/d = n/2) 2HDMs. 

In fig. 1 contours of E^ are plotted. The parameter range -| < yS - i/^^.h < f is sufficient 
as Ed is invariant under the shift yS - i/^j_„ -^ P - ^d,u + ^- As cut-off A = 10'^ GeV has been 
chosen as it is the scenario where maximal FCNCs can be expected and largest deviations of 
the leading log approximation. For this cut-off the grey shaded rectangles are not accessible 



as some Yukawa couplings become non-perturbative below A = 10'^ GeV. As will be shown 
in section 5 E^ can easily be (9(1) while still satisfying the bounds on exotic contributions to 
the FCNCs. 

To evaluate the accuracy of the analytical formulae we have calculated the flavour violat- 
ing Higgs coupling Arf23 solving numerically the full one-loop RGBs in the appendix.- We 
show in fig. 1(b) t^di-hlQdn, and we find a very good agreement with the value Ed calculated 
analytically in eqs. (23) and (25) up to an overall factor of 2.5 (due to the large, flavour inde- 
pendent, eff'ects in the running of the strong coupling constant and the top Yukawa coupling, 
which are not contemplated by the leading log approximation). We find, nevertheless, a new 
feature: There are regions with flipped sign at the top right and bottom left of the figure 
which are shown shaded. In these regions, diff"erences in the running of 7^^ compared to 
y^^\ not present in the leading log approximation, lead to a change of the sign of M^, cf. 
eq. (10). Diagonalizing the quark masses, eqs. (15), this sign is transferred to A^/. 

Before deriving upper bounds on the flavour off"-diagonal couplings we will present our 
approximate formulae in two parameterizations widely used in the literature: the Wolfen- 
stein parameterization of the CKM matrix and the Cheng & Sher parameterization of flavour 
violating couplings in a general 2HDM. 

In the Wolfenstein parameterization Using the Wolfenstein parameterization of the CKM 

matrix [4]: 

( \ - 
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y. 



CKM — 



1 - ^ X AX^ip - IT]) ^ 



-X 1 - f AX^ 



+ 0{t) (26) 



\AX^{\-p-iiq) -AX^ 1 

with X ~ 0.2, A « 0.8, p « 0.2 and 77 « 0.3 [5] and the approximate expressions 

Mf ^- ~ ^diag(/, X\ 1), Mf'- - ^diag(/, /, X\ (27) 
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we get the following estimates: 
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[2^123 ' 


^x\ 


(29) 


and smaller values for the (21), (31), (32) entries. 











We have usedyS - i//e -0 but the result is independent unless cos(J3 - t/Zg) -^ 0. 



In the Cheng & Sher parameterization The non-diagonal couplings of 2HDM are often 
parameterized as [6]: 



Bounds on the coefficients A"', have been derived from experimental results e.g. in [7]. These 
bounds depend on the masses of the Higgs bosons and on whether the parameters A"f are 
assumed to be universal or to posses some kind of hierarchy. We do not assume universality. 
Instead the Yukawa alignment condition leads to: 
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- 6 X 10-^ £„, 


Uial ' 


- 10-4 £„, 


1^31 • 


^ 7 X 10-^ £„, 


(31) 


Ki\ ' 


- 5 X 10-4 Ed, 


I4l ' 


- 6 X 10-2 Ei, 


1^23! ' 


-0.1£,. 


(32) 



5 Experimental Bounds 

In general, there are numerous experimental bounds on the parameters of 2HDMs. Con- 
straints derived for the type I or type II apply also in our scenario as it is more general. In 
[2] bounds on the aligned 2HDM have been studied explicitly. As the alignment condition 
is broken by radiative corrections, in addition tree-level FCNCs are present. This leads to 
further constraints on the parameters of this type of 2HDMs, as FCNCs are known from 
experiment to be highly suppressed. 

5.1 Meson-antimeson mixing 

Stringent experimental bounds on FCNCs come from meson-antimeson mixing. In the SM 
this mixing can occur only at loop level while in a general 2HDM there is also a tree level 
mediation, see fig. 2(a). As the flavour eigenstates are thus not mass eigenstates this mixing 
leads to tiny mass differences that have been determined experimentally for 5°, 5°, D° and 
K^ mesons. Here, we treat only the 5° - 5° system as it gives the strongest constraints (see 
e.g. eqs. (31) and (32)). The eff'ective Hamiltonian of the A5 = 2 transition 5° <-4 5^ is at 
scale ~ mz'. 

Htt^ = Y,CKmz)Ql{mz), (33) 

i,a 

where in a 2HDM with flavour violation at tree level the relevant operators are: 

Ql"^"^ = {bRsdibRsd, Ql"'' = ihsKXhsR), Qf = {bRSL)(bLSR). (34) 



The corresponding Wilson coefficients can be read off" the eff"ective Hamiltonian to be: 



q 



,SLL 



(^^23)' 






m 



A) 



Q - -Kl23^d32 



/ c:2 f2 , ^ 

^ H 



c\ 



SRR 



(A,n2)' 
2 


/ 9 9 


1 \ 

m\ 

A J 


(35) 
(36) 



and the meson-antimeson mass diff"erence can be calculated as: 



Ams, = 



Am™ + \msFl \PfC\\m^) + F^,^^ {c\^\m^) + Cf «(mz))] 



(37) 



where the coefficients P" include both the renormalization group evolution from the high 
scale mz down to low energy ~ mg^ and the hadronization of the quarks to mesons. They can 
be calculated using the formulae in [8] and lattice QCD results from [9]. For the B\ - B\- 
system we get P^ = 3.0 and P^^^ = -1.9. As [Arf],y » [A^]^, for j > i the term involving 
Cl^'^ is always neglible, whereas C^^^ dominates only for degenerate Higgs masses or in the 
decoupling limit {Ca-p -^ 0, mn ~ mA). Therefore: 
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degenerate \ 
masses or 1 . (39) 
decoupling limit/ 



The SM prediction for the mass diff'erence Am^^^ = (135 + 20) • 10"'^ GeV [10] agrees quite 



well with the experimental value Am^ = (1 17.0 + 0.8) • 10^ GeV [5]. Nevertheless there is 
room for new physics as long as the new contribution to the mass diff'erence is smaller than 
the theoretical uncertainty, i.e. 20 • 10^^ GeV. Using eq. (21) this leads to the approximate 
bound (we take F^, = 238.8 + 9.5 MeV [9], me, = 5.37 GeV [5] and values for the quark 
masses at mz from [11]): 
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(40) 



Thus, even for light Higgs masses, (9(100 GeV), the present experimental constraints from 
meson-antimeson mixing allow Eij to be of (9(1). In the case of degenerate masses or the 
decoupling limit the bound is even weaker: 
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(41) 





Figure 2: Feynman diagrams for tree level flavour changing processes in a general 2HDM 

5.2 Leptonic B decays 

The decay B^ -^ l^t (I = e,jx, r), based on the flavour transition b ^> s,h another example 
of a process that can be mediated at tree level in a general 2HDM (see fig. 2(b)) but neither 
in the SM nor in an aligned 2HDM. As the branching ratio depends on the lepton Yukawa 
coupling one expects the decay Bs -^ t^t~ to be favoured. However, the produced r's decay 
immediately to jets and leptons whose observed invariant mass will not reconstruct back to 
the mass of the B meson, so that these decays cannot be tagged in detectors [12]. In contrast, 
tagging B^ -^ iX^jjT is rather easy leading to an unparalleled bound on the branching ratio. 
Litegrating out the Higgs boson, the matrix element for the h exchange is, assuming real Ae22 
(see Feynman rule (74)): 

^5,->/.v- = 4^«-/3(^^23 - A* 32) {B,\ sfb 10) — [—^s„-p + c^-nKii^m (42) 

since (5^1 sb |0) is zero, as B^ is parity-odd, whereas sb is parity-even (see e.g. [13]). On 
the other hand, the matrix element for H exchange is of the same form as eq. (42), with the 
replacements Ca-p -^ Sa-p, Sa-p -^ -Ca-p, mi, -^ mn. Lastly, the invariant amplitude for the 
exchange of a pseudo-scalar A can be inferred from eq. (75) to be: 

^1-.,V- = I^^^'23 + a:,32) <5.vI sy'b 10) ^A,22 /irV- (43) 

The decay rate can now be straightforwardly calculated from the decay amplitudes. Using 

{Brsy'b\0)^ifBms^, (44) 



and neglecting again terms proportional to [/^d]ji,j>i we obtain: 
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(45) 



Requiring that the tree level alone does not exceed the present experimental bound BR(5 , 
fi^fi') < 4.7 • 10-^ [5] we find for A, 22 » v" 
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(600 Gev)2 



(46) 



Hence for Higgs masses of 0(lOO GeV) and E^ of (9(1) there is only a conflict for tanii/r^ - 
/3) ~ l^eii^ <: 50. The opposite case, Ae22 ^ V' results in a bound on E^ and the Higgs 
masses only. However, for generic Higgs masses the bound coming from 5, - 5, mixing is 
stronger. 

6 Conclusions and outlook 

We have studied in this letter the general two Higgs doublet model imposing that the Yukawa 
couplings for each fermion type are aligned at a high energy cut-off" scale. This hypothesis 
guarantees the absence of new sources of flavour violation at tree level. We have shown, 
however, that quantum corrections spoil in general the alignment condition, inducing small 
flavour violating Higgs couplings at low energies. Without imposing ad hoc discrete sym- 
metries and in the absence of tunings, the radiatively induced flavour violation in the Higgs 
sector yields the minimal size of the exotic contributions to the FCNCs in any 2HDM. We 
have shown that, for a wide range of parameters, this exotic contribution is well below the 
experimental constraints from meson- antimeson mixing and leptonic 5-decays. 

In this letter we have concentrated on the constraints in the hadronic sector although the 
analysis can be extended to the leptonic sector. In the minimal extension of the Standard 
Model with one extra Higgs doublet, the Yukawa alignment in the leptonic sector is pre- 
served after including the radiative corrections, due to the exactly preserved family lepton 
numbers. However, this scenario is incompatible with the observed neutrino oscillations. 
In a forthcoming publication [14] we will analyse the implications for the lepton flavour 
violating processes of the 2HDM with heavy right-handed neutrinos, necessary to generate 
neutrino masses through the see-saw mechanism. 
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Appendices 



A RGEs of the Yukawa couplings 

The Renormalization Group Equation of the Yukawa couplings reads 



, dYf 
l6n'fi-^=f3Y,, (47) 



where 



^,,« = a^YlP + J]\3Tr{YlPYir + Yf^Yf)+Tr{Y^^Y^)\Yl^ (48) 



1=1,2 

d ^ d ^ u ' ^ u ^ u ^ u ' 'y d d u r^^u^u ^ u /' 



/=1,2 ^ 



^e,., = UdYf + Y,[^Tx{Y'!''Yf + Y^^Yf)+Tr{Yf-^Yf')]Yf (49) 

'=1,2 

+ Z [-^YfY^!''Yf + YfY^^^^Yf + \Y^^Y^I^^Yf + {Y^^YfY^^^ , 
%, = a.yf + J][3Tr(yry(0 + y<«y(/)t) + Tr(y«y«t)]y« (50) 

/=1,2 
+ ^ /y«y(/)ty(0 1 y(/)y(/)ty(.A 



/=1,2 



11 



where ay (/ = u, d, e) stands for contributions due to gauge interactions, which are flavour- 
diagonal: 

au = -^g]-lg'-^^g'\ (51) 

ad = -8g? - ^/ - ^g'\ (52) 

9 15 

ae = —/ - ^g'\ (53) 

where g.,, g and g' are the gauge couplings constants of SU(3)c, SU(2)/^ and 11(1)7, re- 
spectively. The terms in the first sum in the /3-function are due to the Higgs wave function 
renormalization, the first term in the second sum is due to the vertex renormalization (absent 
for leptons) and the last three (two) terms are due to the renormalization of the fermion wave 
function. 



B Yukawa couplings at the EW scale 

B.l d-quarks 

Yfimz) « 4'V, + ^fYjlY, + Sfy.YlY, (54) 



with 



log — r 
fcj^ = cos i/frf + ^ [ Ud cos ij/d + 'i cos (A* cos((/r„ - il/d)Tv{Yl Yu) (55) 

+ 3 cos ifJd^AYdYl) + cos ifje cos((A: - iArf)Tr(y,y;)], 



log i2z / J 

e^d^ = -TT-T I 2 ^^^ ^d - 2 cos i/r*, cos(iA„ - (A./) | , (56) 

m 31ogf 
^i'^ = ^^cos^,, (57) 



log^r 

kf = sin lArf + -[^[a<i sin i^d + 3 sin i/zl cos((A„ - <Arf)Tr(yJy„) (58) 

+ 3 sin iA./Tr(y^y]) + sin <A, cos(iA: - (Ad)Tr(y,yJ)], 

log ^ / 1 \ 

ef = -j^ I - sin i^d-2 sin i/fl cos(^„ - ^,) 1 , (59) 

n^ 31og^ 
4 =-;:r^sin,A^; (60) 



32;r2 

B.2 u-quarks 



l'®(mz) « -tf^i-^ + ^'^I'rfl'IJ'^ + €^YXYu (61) 
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with 



log — r 

k^J^ = COS (A„ + -r^[au cos (A„ + 3 cos iA„Tr(yJy„) + 3 cos ^d cosii//^ - iAJTr(y^y]) (62) 



1671-2 

+ COS <a: cos((A. - iA«)Tr(y,yJ)], 



log — / 1 \ 

^"^ = ig 2 I 2 *^"^ ^" ~ ^ "^"^ ^'^ cos(i/^,/ - (A„) j , (63) 

^i'^ = ^^cos^„, (64) 

log — 
fc[f ) = sin .A„ + ^^ [ a„ sin (A„ + 3 sin <A„Tr(yJ y„) + 3 sin r, cosiif,, - (AJTr(y,/y]) (65) 



+ 



sin(A:cos(iA.-iA„)Tr(y,yJ)], 



log — / 1 \ 

^u^ = -j^ I - sin i/r„ - 2 sin ip*^ cos((/rrf - (A„) J , (66) 

m 3l0gf 

^L'^ = ^^sin<A„; (67) 

B.3 Leptons 

As there are no vertex corrections in the leptonic sector, the coupling at the electroweak scale 
has a different structure: 

rf (mz) « fcf y. + df^YjX (68) 



with 



log — 
k[}^ = cos (A. + ^^[ «. cos (A. + 3 cos (A: cos((A„ - (A.)Tr(yJy„) (69) 

+ 3 cos lArf cos((A: - iA.)Tr(yrfy]) + cos (A,Tr(y,yJ)], 

31og^ 
5l'^ = ^^cos^„ (70) 



log — r 

kf^ = sin (A. + -j^ [ a, sin i^. + 3 sin i^: cos(iA„ - iA.)Tr(yJ y„) (71) 

+ 3 sin lA./ cos(iA:, - (A.)Tr(yrfy]) + sin iA.Tr(y,yJ)], 

m 31og^ 
5f = ^^sin^,; (72) 
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C Feynman rules 




I 



Mf 



2Cc-p — '-6ij + Sa-p ((A/ + Aj.),-y + (Ay - Ap;yy^) 



V 

for/ = u,d,e 



(73) 




for/ = u,d,e 



(74) 



M/, U/, c/ 





= \ [(^/ - ^P'V + (A/ + A})o-r'] (75) 

for f = d,e and for / = m the same with negative sign 



= -- ^ [(- A^ VcKM + VcKMA,l)ij + (Al VcKM + VcKMA^dijY^] (76) 

for quarks and VcKAiAd ^ A^,, A„ ^ for leptons 

Here m/,- is the corresponding fermion mass and A/ is defined by eqs. (17), (18) and analo- 
gously for the leptons. 
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